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We give a detailed exposition of the "vectorized" notation for dealing with quantum operations. 
This notation is used to highlight the relationships between representations of completely-positive 
dynamics. Vectorization considerably simplifies the analysis of different methods of quantum process 
tomography, and enables us to derive compact representation of the investigated quantum operations 
in terms of the resulting data. 

O I. INTRODUCTION 

(N ' 

Quantum process tomography [l| is one of the standard tools of quantum information science, and efficient methods 
of processing tomographic data are of great practical and conceptual interest. Tomographic data manipulations involve 
massive amounts of the "qubit/qudit algebra" on the finite-dimensional Hilbcrt spaces. Linear operators on them 
qq ■ form vector spaces of their own, and the Hilbert-Schmidt inner product 

A ■ B := tr (A*B) (1) 

i-£h makes such a vector space into a Hilbert space. As a result, all the formulas that describe state evolution and 
'~ H . measurements can be re- written in this vector form. "Vectorized" (or "double ket" Q) formalism is fruitfully used 
' in quantum information (e.g., [H-Q), However, it is sometimes done on the ad hoc basis using not very transparent 
, notation. 

In this article we consistently apply it to various forms of the dynamics of open systems (Sec. Ill) and quantum 
Q" 1 process tomography (Sec. IV). It allows us to give an easy derivation of known results, and provides a number of new 
l—— '■ relationships between the reconstructed data and the investigated processes. To make the discussing self-contained 
| and recognizing that a good notation is a clue to successful derivations we begin with a pedagogical introduction into 
vectorization of finite-dimensional matrices. This is a well-established area of matrix analysis and part of our 
' presentation is devoted to introducing notation which is consistent with other areas of theoretical physics and with 
CO • group representation theory 

in 

II. VECTORIZATION OF MATRICES 
A. Conventions 

o 



We distinguish between upper (contravaraint) and lower (covariant) indices. Components of a vector tf> £ % are 
labeled as ip k , the basis of TL is denoted by {e^}, k = 1, . . . d = dimH. We employ the Einstein summation convention 
5_j ■ where the summation goes over identical upper and lower indices, hence ip = J2k i> ke k = v e k- 

The inner product allows to identify e k = ej . When there is a chance of confusion of the vectorial labels with the 
components, the former are put inside the brackets: the a-th component of the basis vector k is denoted as 



(e fc ) a =e^. (2) 
We choose the orthonormal bases, so — 8^. The Dirac notation is introduced as \k) = e^ and (k\ = e k . Since 

hAH5> fc |fc>, M = 5> k *W> (3) 

k k 

the convention for the components of general vectors is 

W> + ) fe = V fe *- (4) 



Moreover, when the indices are shown explicitly, the sign of a Hermitian conjugation becomes redundant, since 
positioning of the index indicates weather the object is a ket (vector) or a bra (its dual), hence tp k = -0 fc *. 
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Matrices of the size p x q represent the operators between the vector spaces of the corresponding dimensions, 

M:Hi^H 2 , <t) a = M a b i> b , o = l,...p = dim%, 6 = 1, . . . ,g = dim^i. (5) 
If the bases of the two Hilbert spaces are denoted by ei and u& , respectively, then 

m = M^fee' = J2 Mk iW)(i\ g n 2 ®ni (6) 

where is the dual Hilbert space of H. This equation suggests a basis for the space of p x q matrices, 

E k l = \W\, (7) 

which consists of matrices with all but one entry being zero and the unity at the (fcZ)-th entry. The matrix M reads 
then as 

M = M k l E k l , k = l,...p, l=l,...,q. (8) 

Matrices with two lower or two upper indices can be considered as superpositions of direct product states,, as in 

*=F\uieHi®H 2 . (9) 

In writing the matrix element (ab) in any upper or lover index combination, such as M a b , M b or (E k l ) a b , the row 
index a appears first. 

With the above conventions the relationships between components of a transposed matrix are given by 

(M T ) ab = M ba , {M T ) b = M b a , (10) 
the complex conjugation turns the lower into the upper indices, and vice versa, 

(M*) a b = Mt, (11) 

and the Hermitian conjugation satisfies 

(Mt)« 6 = M b : = (M T ) b * = (M*) b a (12) 

Kronicker's delta is real and symmetric, 5 b = 8\ = S b a . Consistency requires to set 5 b * = 5 b . We also keep in mind 
that {E b Y = E b and (E b ) T - (E a b )^ = E b a , so 

M* = J2 M*a b E a b = M a b *E a b , Mt = (M%E a y = ^ M a b *E b a = M* b E b a . (13) 

a.b a.b 

Consider a tensor product of two matrices, K = M ® N. A pair of row indices a\ and a 2 is combined into a single 
row index of K, and a pair of column indices b\ and b 2 is combined into the new column index according to the 
lexicographic ordering rule that is described below. Hence 

n [i — lx 16162 - 1U b! ly 6 2 ' l i4 i 

and similarly for 

K a lJ:=M a b \N b *, (15) 

etc. 

B. Definition and simple properties 

Representation of matrices as vectors on a higher dimensional Hilbert space is called vectorization. It transforms a 
p x q matrix M into a (pq) x 1 column vector denoted by vec(M) or M |. We choose to do it by ordering matrix 
elements lexicographically, i.e. by stacking the rows of M on top of each other to form a vector. Our convention agrees 
with 3], but is opposite to @. We made this choice as to keep the same concatenation rules for both vectorization 
and taking tensor products. 



For example, a 2 x 2 matrix M is re-arranged to form the four-dimensional vector AI , 

^ = {M\,M\,M\,M\) T . (16) 

To automate matrix manipulations we have to lump a pair of matrix indices ab, a = 1, . . .p, b = 1, . . . q into a single 
vector index a. With the stacking convention that we adopted, 

a = f(a,b) = q(a-l) + b. (17) 

The inverse of vectorization restores the matrix form, as mat(A/) = M. The matrix indices are restored according to 

a = \_a/q\ +1, b = amodg = a — q[a/q\, (18) 

where [^J is the largest integer less than or equal to x. 

We now list some useful properties of vectorized matrices [|| . For convenience we consider square matrices that act 
on the space T~L, p = dimH. A matrix 

A = J2 A \\ k )( l \ = A c kc l eH®H*, (19) 

kl 

is transformed to the vector 

A =Y1 Ak i Ck ® Cl = *52( Ac i) ® c * = ( A ® M H c * ® ( 20 ) 

kl I I 

where I p is the identity matrix on T-L. 

Vectorization is obviously linear: for matrices A a and scalars a a , 

vcc(a a A a ) = a a A^. (21) 

Vectorization is intrinsically related to the tensor product. Our stacking convention for a matrix M\ is the same 
as the the index concatenation rule for elements of a tensor product space, such as M ab i— > M a . This will be useful 
in the following. It is also easy to discover how the vectorizations of M T , M*, and M' are related to each other. 
Following the definition of vectorization and Eq we see that 

AT = M\ W = M^+. (22) 



We set M and M* to carry the contravariant indices, so M and M carry the covariant ones. Hence 

M a = M a = AI a b , AT a = M a = M* a 6 = AI a b * = A~I a * , (23) 

and 

AI^ a = AI ta b = AI b *, AI^ a = AI T a b = AI b a = AI^ a *. (24) 
The Hilbert-Schmidt inner product is equivalent to the usual inner product of vectors, 

tr (A*B) = A\ a B b a = = ~A f p ~B*P = = (1, t), (25) 

where j3 = f(b, a). 

Vectorization relates the tensor and the outer product of vectors, 



~2®~3 = ycc{1^), (26) 

as follows from Eq. (f2"2"j) . 

To apply the summation convention we use the opposite label positioning for the basis elements vec(E a b ) = E Q . 
Let a = f(a, b) and /3 = f(k, I). Then the basic definition 

(E> =8i = 5 k a 5 b , (27) 
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and the conjugation rules for vectors and matrices lead, e.g., to 

(^ = (C f )/?=E^=£>\ (28) 

which agrees with 

mat(E*) fc < = {E a b )*t = (E a b ) k ? = {5 k j\)* = 6 a k 6 l b , (29) 

Similarly, 

mat(if )£ = (E a b ) l k = 5 l a 8l (30) 

and 

M = M a if Q , A/t = 1/"*^ = M a if" = if*. (31) 

Consider two matrices A and B, of the sizes p x q and q x r, respectively. Their product C = AB can be written 
as a vector AB* in two ways. The matrix C results from the left action of A on B, so in a vectorized notation 

C a = Ma^B 13 , (32) 

where the matrix Ma will be now determined. Writing the indices in full and using matrix multiplication rules leads 
to 

C\ = M A \ c d B c d = A\B\\\ = A a c I b d B c d , (33) 

so 

M A a b . c d = A\5 b d . (34) 

As a result, 

AB* = (A ® I r )~3. (35) 

Now we see that Eq. (f2T)|) for ~X is a corollary of the above result with B i-> I. The matrix C also results from the 
right action of J3 on A, which is expressed as 

A~B~ = (l p (g> 5 T )^4. (36) 
The following lemma || deals with the triple product of matrices: 
Lemma 1. Let A, B, C and X be p X q, r X s, p X s and (jxr matrices, respectively. Then the matrix equation 

4XB = C (37) 
for gs unknowns X l m is equivalent to the system of ps equations 

(A®B T )1 = (38) 

That is to say, 

Axh = (A®B T )t. (39) 
The proof consists in the repeated application of Eqs. ([35]) and (|36|) . 

Axi = {A® l s )X~i = (A® l s ){l q ® B T )X^ = (A® B T )jt. (40) 

□ 
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C. Reshuffling 

The SWAP operation changes the order of subsystems in a tensor product: 

\i/>)i ® 10)2 SWAP(|V)i ® 10)2) - |0) 2 ® 



(41) 



For two identical systems this operation swaps their quantum states. Consider two vectors i/j € ~Hi and € %2- In 
the component form we write them as 



and 



l,...,p = dim 'Hi, 



(42) 



A single index a of the tensor product \Jj a = (ip < 

a = r(a% - 1) + a 2 . 
Similarly, a single vectorial index for $ = (f> ® xp is 

/3 = q(a 2 - l)m + ai 



a-2 = 1, . . . , r = dimH 2 , (43) 
Q is built from the indices of its subsystems as in Eq. (|18p. 

(44) 



(45) 



The bases of the two tensor products are certain permutations of each other. It is easy to see that this permutation 
j3 = a[r,p\(a) is given by the table 



1 2 

1 1+p 



1 + (s - l)p ... 1 + (r - l)p 



r + 1 r + 2 
2 2+p 



- l)r + 1 



pr 

p+ (r — l)p 



(46) 



where the first row contains indices of the components of ^ that are matched with the corresponding indices of 4> in 
the second row. The vertical lines separate the first and the last r components of VP from the rest, and s = 1, . . . r. 
Its inverse, a[r, p]~ x associates to each element of $ an element of VP and is obtained by interchanging the rows. As a 
result, if we label the SWAP operator as S(r,p), then 



where 



(47) 
(48) 



It follows from this construction that S is an orthogonal matrix, and S(r,p) T — S(r,p) 1 = S(p,r). For example, a 
qubit SWAP is 



5(2,2) 



/ 1 

10 

10 

\0 1 



(49) 



and the SWAP of 2 x 3 systems is done with 



(I 



















1 








1 























1 








1 











\0 











1) 



5(2,3) = 



If we restore the original indices we see that SWAP can be treated as a matrix transposition 

■ bib'2 1 



(50) 



(51) 
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and 

S a2 % 2 =«- ( 52 ) 
For two square matrices M and N that act on the spaces Hi and H2, respectively, the SWAP operation results in 

N ® M = S(r,p)(M ® N)S(r,p) T . (53) 

As we will see in the next sections, it is important to relate vec(M) (g> vec(iV) and vec(M ® N). It is done by the 
operation of reshuffling 0. Let the matrices M and N be of the sizes p x q and r x s, respectively. A matrix element 

M a b \NX (54) 
can be interpreted either according to " first vectorize, then tensor" , as 

M% NX =: C\X ™ ^ ^C A = (vec(M) ® vec(N)) A , (55) 

or as 

M a b \N% =: C a ^ bib2 # C\ n c C B = vec(M ® TV) 5 , (56) 

according to "first tensor, then vectorize" precept. 

A repeated application of the index concatenation definition gives 

A = rsq( ai - 1) + rs(&i - 1) + s(a 2 - 1) + b 2 , B = rsq( ai - 1) + gs(a 2 - 1) + - 1) + 6 2 . (57) 

Hence swapping the indices a 2 and 61 brings M ® iv to M <8 iV. This can be formalized by the following 

Definition 1 (Reshuffling Matrix). For matrices M and TV of the size p x q and r x s, respectively, the reshuffling 
matrix R(M, TV) is defined by 

R(M, /V) ee R(p, g, r, s) = I p ® S(r, g) ® I s , (58) 

Our discussion established the 

Proposition 2 (Reshuffling). Let t/ie matrices M and N be of the size p x q and r X s, respectively. Then 

vec(M ® TV) = R(p, g, r, s)(M ® i^). (59) 

A Mathematica package by Miszczak, Gawron, and Puchala [l(| for the analysis of q uantum states and operations 
implements matrix vectorization and resufling for the use in the relevant functions [l8j |. 

III. APPLICATIONS OF VECTORIZATION TO OPEN QUANTUM SYSTEMS 

The simplest use of the vectorization is in analyzing the convex probability domains of quantum measurements [i[ . 
Consider a positive operator valued measure (POVM) with N outcomes. Vectorization makes it nearly obvious that 
the domain in the probability space formed by probabilities of all outcomes with all possible states of a n dimensional 
system has at most n 2 — 1 dimensions. 

A. Completely Positive Maps 

Evolution of an open quantum system whose initial state is uncorrelatcd with the environment is described by a 
completely positive map. Label the initial state of the system (Hi) as p, and the state of its environment (TL2) by oj. 
If the joint state is given by 712 = p ® oj, and the overall unitary evolution is C/12, then the final state of the system 
is given by 



p ^ T(p) = tru 2 (Ui 2 p®u}Ul 2 ), 



(60) 
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where the operation of partial tracing is defined by 

(tr 2 r)- =r a ^ bia2 . (61) 

Such a map T is trace-preserving, convex-linear and completely positive This latter property means that T 
is positive (maps positive matrices to positive matrices), and also if wc introduce an auxiliary system of arbitrary 
dimension then the map T ® I on the joint system is positive, where I is the identity map on the auxiliary system. 

A map satisfying these three axioms is referred to as a completely positive, trace preserving (CPTP) map. It is 
possible to relax the trace condition to tr [T(p)] < tr [p], allowing for completely positive trace decreasing maps. We 
will only be concerned with CPTP maps and we may refer to maps satisfying all three of these axioms simply as 
completely-positive maps (CP maps). 

Any CPTP map has a convenient operator sum representation: 

Theorem 3 (Kraus representation). A map T acting on density operators of TL is CPTP if and only if there exists 
a set of bounded operators {K n } acting on H such that 

T{p)=^K n pKl where ^K\K n =l. (62) 

n n 

The operators K n are called Kraus matrices and they satisfy K^K n = I, which is known as the completeness 
relation. A Kraus representation of a given process T is not unique. This can be useful as different system-environment 
interactions may still give rise to the same reduced dynamics on the system. 

In general, a linear hermiticity-preserving transformation T (completely positive or not) acting on the space of 
density matrices may be represented by the dynamical (Choi) matrix D(T) pi. [III. Il2j . 

p' a b = D a c /p c d , (63) 

and a vectorized version of this relation 

t a = DV"> (64) 



uses the reshuffled matrix D, 



D\/ = D a c /, (65) 



vec(D) = RD. 

The dynamical matrix has a number of useful properties 0. The trace preserving condition is equivalent to the 
constraint on the partial trace of the dynamical matrix, 

D a c . a d = S d c , (66) 

which implies that the eigenvalues sum up to the system's dimension, ^ Q A a = di. Moreover, if the map is unital, 
i.e., it maps the maximally mixed state into the maximally mixed state, then 

D\ b s =S%. (67) 

The dynamical matrix is Hermitian, D^ a g = D a g, and due to a theorem of Choi [l2[ its positivity is equivalent to 
the complete positivity of T: 

Theorem 4 (Dynamical Matrix). A quantum operation T on a d-dimensional system S is CP if and only if its 
dynamical matrix D(T) is positive- semidefinite (D(T) > 0). 

The proof is based on the eigendecomposition of the Choi matrix, 

D a c b d = D a p = ^2 \nM n a ~Mnl = K~Mn a M* [3 = £ KM^M*^, (68) 

n n n 

where a = f(a, c), f3 = f(b,d), and if all the eigenvalues A„ are non-negative it is possible to define Kraus matrices 
by absorbing the eigenvalues, 



Kn = VKM n . (69) 
We also get a compact expression for the reshuffled matrix, 

D = (70) 
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B. Linear Superoperator 

Vectorization allows a different perspective on Choi matrix and its reshuffled version. We consider a combination 
of basis matrices E b and T(E a b ) = E' b , with suitably arranged component indices. 
A linear superoperator $j acts on vectorized density matrices and is defined by 



On the one hand, since 



we have 



As a result, 



$ T t = T(p). (71) 

<E Q , E^> = E^eJ = tr (E a ^E c d ) = S a c S d = 8% (72) 

p^E^E^E^ = /E^ = $ T ^. (73) 

$ T =Y,^L^ =K^- (74) 



A comparison with Eq. (1641) identifies the linear superoperator with the reshuffled dynamical matrix, $j = D(T). The 
matrix elements of $t are obtained by 

^ = e1^ x = J2(E' a b ) k i(E a b *)£ = {E' a b )\5 a m 5Z = {E'-)\ = {E'^y. (75) 

a . b 

Indeed, p f K = so 

0\ m n = (E'-)\ = J2(K x ) k m (K* x )r. (76) 

X 

We note that when expressed in terms of the Kraus matrices K n the reshuffled matrix D is built according to "first 
tensor, then vectorize" prescription. However, the opposite approach "first vectorize, then tensor" is taken when $j 
is expressed in terms of the basis matrices and their transforms. 

C. Process matrix 

We define the process matrix Aj for a map T as a matrix on a fictitious space Hi <g) %2, where the two spaces are 
copies of H. Without going through the vectorization of Eqs. (|72[) . ([73")l . the transformation p' = T(p) can be written 
as 

p' = J2 P C d E'a tr (E a b ^E c d ) = p c d E' a b tr (E b a E c d ). (77) 

a . b 

Thus it is given by 

T(p) = tr 2 [A^ 2 (I ® p)] , Ay 2 ee E' b ® E b a . (78) 

For our purposes it is more convenient to use the matrix arrangement of $t, and define a new four- indexed object 
according to the "first tensor, then vectorize" prescription, 

A T ee E' b ® E^ = £ E' a b ® E b \ (79) 

a . 6 

where we keep a complex conjugation as a reminder for the correct index placement. Then 

T(p)=tr 2 [A T (I®p T )], (80) 

and 

At = (£ E' b ® E a b *)\ m n = (E' a b ) k m 8?6Z = (E'^m- (81) 

a,b 

A comparison with Eq ((55)) shows that D(T) = Aj. 
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D. Summary of relationships 



The Jamiolkowski isomorphism 0, [l3| identifies dynamical matrices of CP maps with certain entangled states. 

Consider a generalization of the Bell state |$ + ) = (| + +) + | ))/2 on the space T-Li ® H2 of the previous section. 

Its density matrix is given by 



r + = J2E a b ®E a b /d 1 



(82) 



where d = dim'H. It is now easy to write the matrix elements of p' = (T ® \)p, 

r' + k \ mn = (E'n k Jd = A T *U"M 

which establishes 



(83) 



Theorem 5 (Jamiolkowski isomorphism). Any linear map T acting on the space of mixed states on the Hilbert space 
% can be associated, via its dynamical matrix D(T), with an operator on the enlarged Hilbert spacee 'H^'H, 



d(t) = (t®i)t+ = t;, 

where for a completely positive map T the matrix t' + is a valid quantum state. 



(84) 



As a computational aid we present a compact summary of the relationships between different representations of a 
completely positive map T. 



To\From 


D(T) 




K n 


T(p) = 


tra [D(Ii®p T )] 


mat($T~p*) 


^ n KnpKt 


D = 


D 


mat(R _1 $^) 




$T = 


mat(R~D) 






Kn = 


\/\ n ~mat(Mn) 







TABLE I: Relationships between mathematical representations of CP maps. 

Two possible viewpoints on both D and <E>t as obtained either according to "first tensor, then vectorize" or "first 
vectorize, then tensor" approaches (with a counterpart being obtained by the other procedure) are consistent, because 
the reshuffling matrix involved in the process satisfies R = R _1 . The latter result follows from Eq. (|46|) . with all the 
dimensions being equal to dimH. 

IV. PROCESS TOMOGRAPHY USING THE SUPEROPERATOR REPRESENTATION 

A. Dual bases 

Quantum process tomography is used to characterize the properties of quantum information-processing devices. 
The standard process tomography [l[ uses a tomographically complete set of input states that spans the entire set of 
states. The input states are sent through the investigated device and its action is then reconstructed by analyzing 



the outputs. A set of dual states D = {D^} for a tomographically complete set p m = {p^}, p = 1, 



d 2 , where 



d = dimH, plays an important role in the mathematics of reconstruction. It is defined by the orthogonality relation 
ti (Dlp^ = 8 Vfl , V[i,v. 

If we rewrite this in a vectorized notation we obtain 



(85) 



hence 



(86) 
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Introducing the operator 

P^P^^Ip 1 , (87) 



we see that PD V = p v , hence 



It is possible to write this relationship between a tomographic set and its dual in a simpler and computationally 
more efficient way. We write vectors of the set p m as the columns of a matrix 

[Pin] = {pi,...,pt)- (89) 

Since its columns are linearly independent, this d 2 x d 2 matrix is invertible. Its Hermitian conjugate can be written 
as 

; 

With this notation we get 



ft 



^=[Pin][p i „] t . (91) 

As a result, 

K =[p in ]- lt [p in ]- 1 ^. (92) 
If we introduce [D] = (-Di, . . . , D n ) } then an even simpler expression is obtained: 

[D} = [pi n ]-^[p- m ]- 1 [p i n} = [ Pi n]-^. (93) 

Now we apply this to the standard process tomography. 

B. Process tomography with $y 

Proposition 6. Let p; n = {p^} be a tomographically complete set of input states with dual basis D = {D^}. Then 
the linear superoperator $y for a CP map T : p i— > T(p) is given by 

— > 



^t^TW^^;^ (94) 

Proof. Consider an arbitrary state p = p^p^. We have to show <&j~ff = vecT(p). Applying the above expression to 
p we see that 

*t(7) = P"Z (D„ pv) = p v ZW = P^Z = T(p). (95) 

□ 

We simplify this expression by using the matrix of vectorized input states, Eq. (|89j) . This leads us to the following 



Proposition 7. For the set of output states p Q ut = {T(/9 M )}, where p ln = {p^.} is a tomographically complete set of 
input states, the linear superoperator <&j is given by 

$T = [PoutHPmr 1 (96) 

Proof. We use Eq. (j94|) . 

[Poutlbin]- 1 = kut]^ =^ct'=^j. (97) 

□ 
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We rewrite this result in terms of the probabilities of various experimental outcomes. 

Proposition 8. Let {M^} be a tome/graphically complete measurement set (M^ > 0, J~] =1) with a dual basis 
{E v }, tr(EjAf^) = <^. Then 

$ T = [E] [m] [£] f , (98) 

where [E] = (E\, . . . , E n ) is a matrix of the vectorized dual elements, and [to] = (m^) is a matrix of probabilities, 
m ia ,=tr(M*f/ v ). 

Proof. Since any p is reconstructed according to p = E^m^, where to m = tr (M^p), 

[Pout] = [E][m], (99) 
and [D]t = [pij- 1 . □ 

C. Ancilla-assisted process tomography with linear superoperator 

Presentation of Jamiolkowski isomorphism and the manipulation of data in the ancilla-assisted quantum tomography 
(AAPT) @, M, 

also benifit from the vectorized notation. 
We introduce an auxiliary system (ancilla) H2 to our principal system "Hi, so that the state space of the joint 
system is given by Hi ® %2- AAPT aims to reconstruct the CP map T on the states of T-L the action of T ® I on a 
single state T12 of this combined system. Evolution on Hi results in an operation T described by a superoperator $j, 
and the ancilla does not evolve. Any initial state T12 of the joint system can be represented as 

712 = $^ Y. W>1 = 1 ( 10Q ) 

which is in general entangled (it is separable if and only if all u/ 1 > 0). 
Lemma 9. The joint system dynamical matrix has the form 

$t®i = R($T®T)Fr\ (101) 

where R is the reshuffling matrix, which in this case again satisfies R = R -1 , and the linear superoperator of the 
identity map is I = $|. 

Proof. Linearity of the evolution is expressed as 

r' = (T0l)r = ^Kj"Ty«w M . (102) 

Hence 



=X) W V,.®^ = ^^R(^'«>^) = R($t®T)^^^®c^. (103) 

MM M 

Matrix elements of T satisfy T fe z m n = 5^ n 5\ l = I a p. Another reshuffling leads to the desired result, 

"t*' = R($ T «)T)R- 1 ^. (104) 

□ 

Both vectors 

^ A = = (R- 1 ^)^ = (R- 1 ^ m n = w*pA <V£, (105) 

M 

and t ou I a = (R~ 1 7'') Q ' 9 correspond to the right hand side of Eq. ([5^1 . Their relationship is of the forms of Eq. (|3"5j) : 

= ($ T «^. (106) 

Introducing $ T = mat(ri n ) we have 

*^ = (*T®I)7£ = (I® (107) 

Hence we can recover the dynamical matrix $y from the output state tas when the matrix $ r is invertiblc: 
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Proposition 10. A linear superoperator $y is recovered from the output of AAPT with the initial state r according 
to 

^ = (I®(<&-T)Vy^ 7<5 ( 108 ) 
An important special case is the entanglement-assisted process tomography, where the input state is a maximally 
entangled 

T = T+ = j2\m®\m/d. (109) 

It corresponds to $ r+ = mat(7+) = I. Hence we established a useful expression for the dynamical matrix and a dual 
form of the Jamiolkowski isomorphism 

Corollary 11 (entanglement-assisted process tomography). In the entanglement assisted process tomography (AAPT 
with the maximally entangled initial state t+ ) the dynamical matrix the linear superoperator "fy is determined by the 
output state according to 




V. CONCLUSIONS AND OUTLOOK 



Adopting vectorized notation allows transparent and consistent representation of various forms of open system 
dynamics, isomorphism between states and operations and representation of process tomography. Neat expressions 
for process tomography ([551) . (|108|) and (|110p use already reconstructed output states. Processing of actual state 
to mog raphic data is much more involved. In particular, relative frequencies cannot be directl y ta ken as probabilities 
[J, Il4| . assumption of completely positive dynamics should be justified or may not be true [l5l.ll6j , matrices may have 
only generalized inverses [17] . A transparent and versatile notation is a great asset in dealing with these issues, and 
we expect that it will simplify some of the existing formulas and bring to light new useful relationship, similarly to 
presented in this work. 
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